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ABSTRACT 



This thesis supports computer-aided data analysis of holograms produced from 
rocket motor firings. The work reported in this thesis modified existing software code to 
make it compatible with installed upgrades in the microcomputer imaging system. In 
particular, this involved converting the format of C language function calls to ITEX/PC 
image processing software to that dictated by ITEX/PCp/us software. Additional modi- 
fications were performed to enhance code portability and optimization. Results indicate 
that all efforts to incorporate the new system software have been successful. 
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I. INTRODUCTION 



A. BACKGROUND 

This thesis was a continuation of the work currently being performed at the 
Naval Postgraduate School and was supported by the contributions attained in Ref- 
erences 1, 2, 3 and 4. The goal of this research was to obtain data concerning the 
effect of propellant properties upon the performance of solid propellant rocket mo- 
tors. Data analysis is performed on holographic images obtained from combustion 
products created during rocket firings. Holograms are obtained through the use of 
a pulsed ruby laser and a holocamera. Reconstruction of the holographic image for 
analysis is achieved with a krypton laser, variable-power microscope, and video cam- 
era. Unfortunately, the optical diffuser required to record the hologram creates a 
speckle-corrupted image which must be filtered before the data retrieval process is 
initiated. The data attained from these procedures is used in performance studies, 
stability analysis predictions, and exhaust signature diagnosis. 

B. THESIS OBJECTIVES 

The thesis objectives were: 

1. To implement software modifications in existing local code in support of the 
new computer system configuration. 

2. To enhance program code portability. 

3. To optimize program code and attempt to reduce user interaction where possi- 
ble. 

4. To create a menu format for executable files. 

All objectives were accomplished. ITEX/PCp/us library routines have been suc- 
cessfully integrated into the image analysis programs. The program file genfunc.c was 
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used to optimize program size, and new functions were created for user convenience. 
A menu format was devised to make the image analysis procedure easier to use. 

C. SOFTWARE SUPPORT 

The associated system software consists of the following items: 

• Image analysis programs created via thesis research in References 1, 2, 3, and 

4, 

• Imaging Technology’s ImageActionp/us Software package, 

• Imaging Technology’s ITEX/PCp/us software library, 

• IBM Disk Operating System version 3.3 (DOS), 

• Microsoft C optimizing compiler version 5.0 (MSC), and 

• Microsoft Codeview version 2.10. 

The ImageActionp/ws software package is a menu-driven system that allows the 
user to process stored or live video images [Ref. 5]. The two- level menu organization 
provides image analysis, graphics, geometric operations, and image processing capa- 
bilities. This software package also has the ability to utilize user-generated special 
purpose script files. The script files are created in C language and are supported by 
libraries created by Imaging Technology, Inc. Imaging Technology’s ITEX/PCp/us 
software libraries are compatible only with the Microsoft PASCAL, FORTRAN, and 
C compilers. These libraries consist of image processing and graphics function sub- 
routines used in conjunction with the PCVISIONp/us Frame Grabber board, a frame 
memory and video digitizer installed in the IBM PC/AT. 

D. THESIS ORGANIZATION 

This thesis is organized into five chapters. The introductory chapter provides 
the project background and goals. It also discusses thesis objectives and organization. 
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Chapter II describes the system hardware modifications. The PCVISION and 
the PCVISIONp/us hardware configurations are compared, and the respective hard- 
ware capabilities are presented. Also the hardware performance is briefly discussed. 

The third chapter gives a synopsis of the image processing procedure and dis- 
cusses actual software programs and improvements at length. The chapter also de- 
scribes the techniques and software development tools utilized for program modifica- 
tion. Finally, Chapter III provides program specifics. 

The software performance analysis is provided in Chapter IV. Program analysis 
consisted of comparisons in the following areas: executable program size, program 
execution times, and program result verification. 

Lastly, Chapter V summarizes the findings obtained during the course of this 
research and provides recommendations for future thesis research. 
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II. SYSTEM HARDWARE MODIFICATIONS 



A. HARDWARE CONFIGURATION 

Digital image processing involves the achievement of the following objectives 
through a computer interface: recognition, segmentation, enhancement, or analysis 
of an image. The image processing and digitization system consists of the following 
components: 

• An IBM PC/AT computer with 40 megabyte hard drive, 

• An Intel INBOARD 386/AT microprocessor and 80387 math coprocessor, 

• An AST ADVANTAGE memory expansion board, 

• A PCVISIONp/us Frame Grabber board, 

• A computer monitor with EGA display, 

• A Panasonic color video data monitor, and 

• A video cassette recorder. 

A block diagram of the basic elements of the image processing system is presented 
in Figure 2.1. One of the primary components of this image processing system is 
the video cassette recorder, which provides the introduction of the image into the 
PCVISIONp/us Frame Grabber board. This image signal consists of timing and 
analog video information in a standard RSI 70 video format. Through the digitization 
process, the analog signal is converted to a digital format for processing and storage. 
The IBM PC/AT central processing unit accesses these digital values via the frame 
memory. 

The PCVISION plus Frame Grabber is a single board that is connected directly 
into one of the 16-bit expansion slots of the IBM PC/AT computer. Pixels are stored 
in frame memory after the Frame Grabber digitizes the analog video signal at a 
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Figure 2.1: Elements of the Image Processing System. The Dotted Rect- 
angle Encloses Frame Grabber Board. (From Ref. 7) 
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Figure 2.2: Display of AVW and Hidden Areas of Dual Memory Configu- 
ration. (From Ref. 6) 

rate of thirty frames per second with eight hits of accuracy. In previous work, a 
PCVISION Frame Grabber board was utilized (versus the PCVISIONp/us version). 
The PCVISIONp/us board permits either two 512 x 512 pixel images or one 610 x 
512 pixel image to be stored within frame memory [Ref. 7:p. 1-4]. In the PCVISION 
board, only one 512 x 512 pixel image can be stored [Ref. 8:p. 1-2]. The bottom 
thirty-two lines of the video output from both versions of the Frame Grabber board 
are within a “hidden area” which can be accessed by the host computer but not viewed 
on the video monitor without scrolling the image. The active video window (AVW) 
is the area within frame memory in which the pixel values are displayed or stored via 
the scanning circuit [Ref. 6: p. 1-14]. Figure 2.2 and Figure 2.3 indicate respectively 
the dual and single memory configuration hidden and AVW areas. Currently, only 
the dual memory configuration is utilized in thesis research. 
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Figure 2.3: Display of AVW and Hidden Areas of Single Memory Config- 
uration. (From Ref. 6) 

Prior to display of the image on the video monitor, look-up tables (LUT) are 
employed for the transformation of pixel intensities. The original PCVISION Frame 
Grabber had only four LUTs associated with output [Ref. 8:p. 2-8]. The current 
Frame Grabber board has one input digitization path and three output channels, each 
allocated eight LUTs [Ref. 7:p. 1-4]. These LUTs allow for real time processing by 
virtue of their ability to afTord simple point transformations of the 256 gray levels 
without any processing delay. Look up tables are addressed by eight bits of pixel 
information provided by the frame memory. A digital- to-analog converter and the 
look up tables constitute the display logic unit of both versions of the Frame Grabber 
boards. The digital-to-analog converter receives gray intensity values corresponding 
to the actual value addressed within the LUTs [Ref. 8:p. 2-9]. 
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Both versions of the Frame Grabber board are compatible with the IBM PC/AT 
computer bus structure. This host interface provides the communication path be- 
tween the Frame Grabber boards and the personal computer. The previous PC- 
VISION board host interface consisted of three components: the Interrupt Logic, 
Control Register Interface, and the Frame Memory Interface [Ref. 8:p. 2-1]. The 
PCVISIONp/us host interface is composed of the Control Register Interface and the 
Frame Memory Interface [Ref. 7:p. 4-2]. The Interrupt Logic had afforded the 
PCVISION Frame Grabber the ability to be controlled in real time while the IBM 
personal computer processed other instructions [Ref. 8:p. 2-3]. This outdated design 
feature has been superceded by transparent computer access to frame memory in the 
PCVISIONp/us system [Ref. 7:p. 1-6]. 

Within the input/output (I/O) space of the PC/AT, sixteen bytes are reserved 
for the PCVISIONp/us Frame Grabber control registers. Currently, only twelve of 
these sixteen bytes are being used by the system. Access to the remaining bytes 
is prohibited; they are reserved for future system growth [Ref. 7:p. 4-3]. On the 
contrary, the PCVISION board requires 32 bytes within the I/O space for its control 
registers. Only the first seven of these bytes are actually utilized for control functions 
[Ref. S:p. 2-3]. 

The frame memory of the PCVISIONp/us system is apportioned into eight equiv- 
alent blocks of 64I\ bytes. Block dimensions are 512 pixels by 128 pixels, and the 
physical layout of the dual-store configuration frame memory is showm in Figure 2.4. 
This configuration can be contrasted to that of the 640x512 image in which memory 
is apportioned into eight blocks of 1024 pixels by 64 pixels. The single-store frame 
memory configuration can be viewed in Figure 2.5. Individual selection of these blocks 
is controlled by select bits within the Control Register. As viewed by the host, each 
block gives the semblance consecutive data bytes originating at the memory base 
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Figure 2.4: PCVISIONp/us Dual Frame Memory Layout. (From Ref. 7) 



address [Ref. 7:p.4-3] . For the PCVISION system, the frame memory was subdivided 
into four blocks of 256 pixels by 256 pixels, forming a quadrant structure as displayed 
in Figure 2.6. At any time only one block can be accessed from the personal computer 
bus [Ref. 8:p. 2-7]. 

Both versions of the Frame Grabber board require the assignment of a single 
block of 64K bytes within the IBM PC/AT memory space located above 640K. Spe- 
cialized software controls the mapping of individual blocks of frame memory into this 
reserved space [Refs. 7 and 8]. In the case of the PCVISIONph/s board, the memory 
base address is factory-configured at address A0000II. This address may be modified 
upon board installation provided the subject address selected is within the guidelines 
set forth in Reference 7. 

It was originally intended that this thesis work would be implemented on a 
COMPAQ DESKPRO 386/20 computer system. The PCVISIONph/s memory base 
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Figure 2.5: PCVISIONp/us Single Frame Memory Layout. (From Ref. 7) 
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Figure 2.6: PCVISION Frame Memory Layout. (From Ref. 8) 
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address was configured to an approved value of D0000H. Therefore, the address region 
from D0000H to DFFFFH was reserved within the computer memory space, and all 
other memory and peripherals were restricted access to that memory region. Upon 
subsequent testing and diagnostic checks, it was found that there was an unresolved 
memory conflict within the COMPAQ computer memory which affected some of the 
operations of the Frame Grabber board. Due to necessity, the board was consequently 
installed in the IBM PC/AT. Memory base address D0000H was again selected; this 
time there were no resident memory conflicts. 

B. HARDWARE PERFORMANCE 

With the introduction of the PCVISIONp/us Frame Grabber board, the ability 
to process two separate images in frame memory is invoked. The dual memory con- 
figuration allows the host to select either image area for analysis and manipulation. 
A display memory function allows the user to view independently either image area. 
These capabilities allow for considerable future flexibility in the analysis of the rocket 
motor holograms. 

The PCVISIONp/us board not only provides more image processing capabilities 
than the PCVISION board, but it also establishes the need for additional function 
libraries to support its abilities. For this reason, image analysis programs devised in 
prior thesis research had to be modified to reflect the new supporting libraries. These 
programs and their resulting modifications are emphasized in the following chapter. 
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III. SOFTWARE PROGRAMS AND 
IMPROVEMENTS 



A. IMAGE PROCESSING PROCEDURE 

As was discussed in the introductory chapter, holographic techniques are em- 
ployed to obtain the image of combustion products from small rocket motors during 
firing. During the hologram reconstruction process, a krypton laser is used in con- 
junction with a variable power microscope. Due to the presence of a diffuser that is 
required to eliminate thermal gradient effects, speckle noise is introduced into the cap- 
tured image. The following paragraphs describe the procedure developed to acquire 
and process this image for statistical study. 

1. Image Acquisition 

The speckle-corrupted image is recorded for future processing via a combi- 
nation of a video cassette recorder and low-light-level camera attached to the variable 
power microscope. The magnification power of the microscope determines the rel- 
ative size of the pixel elements comprising the image. Pixels are the standard unit 
for the measurement of the features within the image. Object sizes are determined 
by scaling the measured object width (in pixels) by the dimension spanned by one 
pixel (in microns per pixel). A quantization error is introduced into this calculation 
due to the necessity of measuring dimensions with an integer number of pixels. As 
the width decreases, the proportional level of uncertainty increases, which produces 
a more erroneous scale factor calculation. A calibration object is inserted within the 
holographic image to determine the width spanned by a pixel. It is essential that this 
object covers an appreciable portion of the monitor screen (at least one half of the 
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area) to counteract this quantization error. The calculation of the conversion factor 
from pixels to microns becomes increasingly important in subsequent steps of image 
processing. 

2. Image Digitization 

From the VCR tape obtained from the hologram reconstruction, the menu- 
driven ImageActionp/us software is used to “grab” and digitize the desired image. 
The PCVISIONp/us Frame Grabber board implements this procedure. The digitized 
video image is comprised of a 512 x 480 array of pixels having gray intensity levels 
ranging from 0 to 255. These intensity levels extend from blackest-black, level 0, to 
whitest-white, level 255. 

3. Speckle Reduction Filtering 

Image filtering is the paramount technique used in speckle reduction. The 
amount of speckle reduction achieved by filtering algorithms is gauged by a calculation 
of the speckle index. Speckle index is defined as the ratio of the image’s local deviation 
to the local mean of pixel values, and this measure is indicative of the random speckle 
noise inherent in the image [Ref. l:pp. 17-19]. 

The three filtering algorithms are based upon synthetic aperture radar 
(SAR) theory presented in References 9 and 10. There are three separate filter mod- 
els. Derived from nonlinear techniques, two filter designs are based upon statistical 
methods. The third filter model has its foundation in a geometric hulling algorithm. 
The geometric hulling algorithm produces the best results by enabling the user to 
discontinue filtering at a less severe level of resolution degradation. 

4. Image Threshold 

In an ideal image, feature particles and the background would have distinct 
gray intensity levels. The process of thresholding creates a binary image in which 
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image features will appear black and the background white. In reality, the thresholded 
image will contain some amount of speckle manifested as feature particles and will 
have lost portions of actual particles during the filtering and thresholding process. 
Image evaluation errors are inherent whenever the grey levels of the object overlap 
the gray levels of the image background. The selection of a threshold value is obtained 
from observation of the filtered image. ITEX/PCp/us software allows the user to select 
a threshold value for which the designated LUT is modified. This modification is not 
made permanent until the image is mapped into frame memory, allowing the user 
to rapidly iterate values to make a determination on the most appropriate threshold 
value. 

5. Feature Identification 

Once the filtered image has been thresholded, it is ready to undergo fea- 
ture identification. It is imperative that the programmed image has been previously 
thresholded due to the logic of the identification algorithm. This procedure locates 
all connected adjacent feature pixels and assigns feature value numbers to them. 
Through this labeling process, interconnected pixels are identified as an individual 
object. The total number of objects recognized serves as input for the next step 
of image processing. This number is restricted to 32,767 or below based upon the 
limitations imposed by the C language on integer values [Ref. 4:p. 9]. 

6. Feature Sizing 

The procedure of feature sizing determines both the maximum vertical and 
horizontal diameters of the identified objects. The total area of the subject particles 
is computed, and all these measurements are delineated in terms of microns. As 
described before, a conversion factor is necessary for manipulation of the feature 
particle’s pixel measurements into actual physical dimensions. This constant, called 
SCALE-FACTOR, is presently based upon the dimensions of a threaded screw at the 
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same magnification of the captured image. The screw threads serve as a calibration 
object to transform image object pixel size data into actual physical measurement 
criteria. The final output of the feature sizing routine is presented in tabular form 
suitable for input into the STATGRAPHICS software package, a statistical analysis 
tool. 

7. Histogram production 

To exhibit feature size data in a histogram form for ready evaluation, the 
STATGRAPHICS software package is utilized. This particular software package was 
selected because it is compatible with the IBM-PC/AT system and is relatively easy 
to use. Histograms produced by STATGRAPHICS are used for further statistical 
analysis to determine total particle distributions of holograms. 

B. PROGRAM OPTIMIZATION 

The primary purpose of this thesis research was to modify existing C language 
routines to make them compatible with the ITEX/PCp/us software. Subsequent pro- 
gram optimization was performed in support of this goal. One significant improve- 
ment was the removal of all general-use functions from the main programs. Although 
a general support functions program file had been created in the previous rendition, 
this file had not been adequately utilized. As a result, all the main programs had 
needlessly large file sizes. 

The C language provides the ability to use the infinitely-abusable GO TO state- 
ment. This creates a situation in which branching occurs in a deeply nested structure. 
The GO TO organization may prove effective in some instances, but it always pro- 
duces confusion for the individual having the misfortune to try to modify that code 
sequence. For this reason, program structures were modified to eliminate all GO TO 
branch statements, thereby making code more modular and portable. 
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A menu format was developed for use with the executable programs. This 
format will facilitate demonstrations of image processing and makes the routines 
“user- friendly”. In addition, supporting functions were devised for some of the pro- 
grams to further eliminate the need for user interaction with the programs once the 
selected image had been loaded and image processing regions were identified. These 
modifications facilitated more rapid program execution. A more detailed discussion 
of program modification specifics will occur in a succeeding section. 

C. SOFTWARE DEVELOPMENT TOOLS 

The Microsoft CodeView Debugger program was utilized for program develop- 
ment. It requires that specific compiler directives be used during program compi- 
lation. Although Codeview is an excellent developmental tool, program debugging 
can become a tedious process. No program modifications can be performed while in 
Codeview, and for that reason, it is not an adequate tool for what is needed for pro- 
gram modules of this size. The ITEX/PC plus routines require large memory module 
handlers, and few other debuggers are able to meet this requirement. Quick C version 
2.0 is capable of handling large memory modules and has an online source debugger 
which is more versatile than Codeview [Ref. 11]. It is recommended that in future 
work with these program modules, Quick C should be included in system upgrades. 

Several batch files were created for compilation and linker processes since Code- 
view required special compiler directives. There are two different versions of these 
files: compilation and link with debugger information, and compilation without de- 
bugger information for file optimization and compaction. Listings and justifications 
for these DOS Batch files are referenced in Appendix K. 
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D. PROGRAM SPECIFICS 



The succeeding paragraphs discuss the programs that were developed in support 
of this thesis research. These existing programs were modified to be compatible with 
the ITEX/PCp/ws software libraries and the ImageActionp/us routines. Where at 
all possible, program optimization was performed. In some cases, code had to be 
completely revamped to accommodate new requirements. All source code listings are 
included as appendices at the end of this document. A block flow diagram of program 
usage is provided in Figure 3.1. 

1. Program File thesis. h 

The program header file thesis. h (Appendix A) contains all requisite library 
include files, C manifest constants, and MACRO definitions as necessitated by main 
programs. For code portability, it is imperative that a centralized copy of all common 
declarations and definitions be maintained as the programs evolve. Therefore, any 
modifications to common parameters can be performed by substitutions within this 
file. 

With the advent of ITEX/PCp/ws software, additional include files were 
necessary beyond the scope of the original header file listed in Reference 4. For 
example, ITEX/PC libraries require the use of the include file itexpc.h [Ref. 12:p. 
1-15]. This requirement has been superceded by the ITEX/PCp/ws software which 
dictates the use of itexpfg.h and stdtyp.h include files. Not only must these include 
files be listed in the header file, but to be utilized, they must also be resident in the 
actual program directory [Ref. 6:p. 1-9]. 

The added capability of dual memory required that many of the calls to the 
software libraries be altered. Additional manifest constants were needed for memory 
specification. They included MEMORY, defining the type of frame memory con- 
figuration, and constants necessary to define the size of the selected frame memory 
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Figure 3.1: Flow Diagram of Image Processing Programs. 
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configuration: XSIZE, YSIZE, and DEPTH. XSIZE and YSIZE indicate respectively 
the horizontal and vertical coordinates of the memory size, and DEPTH refers to the 
number of bits assigned per pixel in frame memory [Ref. 6:p. 2-6]. 

2. Program File geofil.c 

Appendix B contains the listing for geofil.c , a program based upon a geo- 
metric hulling algorithm introduced in References 9 and 10. The algorithm is applied 
to the image and its complement by performing comparisons between a central pixel 
value and adjacent pixels in horizontal, vertical and diagonal directions. The original 
geometric filter routines were adapted by Edwards [Ref. 1] from information set forth 
in References 9 and 10; and Kaeser successfully adapted these to the C language [Ref. 
4]. Edwards’ FORTRAN filtering routine was capable of handling only a quarter of 
the actual image due to memory restrictions, whereas the C language version fully 
supported the entire image. A large part of the C version’s success was due to the 
utilization of virtual memory arrays for data processing. These arrays are stored on 
disk, but they are accessed as though they are stored in memory. File constraints 
within the associated operating system po^e the only limitation upon the virtual array 
size [Ref. 13]. 

Few modifications were performed on the geofil.c program to make it com- 
patible with the ITEX/P Cplus software. The majority of the changes for compatibil- 
ity encompassed general support functions. However, much of the C language code 
was altered to attain optimum modularity. This filter program had incorporated the 
use of GO TO statements to facilitate branching between the filter stepping functions. 
Alterations included the removal of these subject statements and the formation of two 
recursive functions. These recursive functions support program modularity and per- 
mit easier modifications for future updates. 
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One common modification performed on geofil.c and its other main pro- 
gram counterparts involved the declaration of the image filename and comment line 
(comline). The header file used by ITEX/PC, itexpc.h, required special file type struc- 
tures used in conjunction with parameter passing for SAVIM and READIM functions. 
These structures consisted of strings with length 20 or 200 characters, and were clas- 
sified as typedef LS20 or LS200, respectively [Ref. 12:p. 1-19]. ITEX/PCp/us limited 
filename to 20 characters and comment lines to 200 characters, but removed the 
special structure requirement. SAVIM and READIM routines now require that the 
image filename and comment line be classified as pointers of type character [Ref. 6]. 
Consequently, all program files were altered to reflect this requirement. 

3. Program File Istat.c 

The program Istat.c (Appendix C) uses local statistical methods to calcu- 
late a weight, k, which determines the gray level of a pixel [Refs. 9 and 10]. The 
mean and variance local estimates within a 5x5 window encompassing the selected 
pixel determine the value of k. The k value then specifies the new gray level of the 
pixel in question [Ref. l:p. 24]. Virtual arrays again assist in the data manipulation 
and storage processes. 

Kaeser's version of this filtering algorithm required that the user input the 
standard deviation of the image directly into the program [Ref. 4]. The operator had 
to obtain this value from the histogram function of the ImageAction software [Ref. 
14:p. 11-9-3] before running the filter program. ImageActionp/us software no longer 
provides the image mean and standard deviation calculations [Ref. 5]. Consequently, 
a new function, dev(stddev). was devised to calculate the image standard deviation 
and input it directly into the filter routine without user intervention. This new' 
function presents a vast improvement in terms of user convenience. 
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4. Program File 2sigma.c 

The 2sigma filter algorithm provides the basis for 2sigma.c (Appendix D) 
and was adapted from the program introduced by Edwards [Ref. 1]. From the 
standard deviation of a Gaussian distribution, it is known that 95.5% of a selected 
region’s pixels will be within two standard deviations of the local mean. An averaging 
formula has been devised in which pixels meeting the “2sigma” criteria are averaged 
for a given pixel, and those outside the range are rejected [Ref. 10]. These calculations 
are performed for 5x5 arrays of pixels. The program 2sigma.c utilizes virtual arrays 
and the standard deviation function as well. It also has reduced size due to the 
removal of all general supporting functions into the program file genfunc.c. 

5. Program File threshit.c 

Threshit.c (Appendix E) produces a binary version of the filtered image. 
The program operator enters a threshold value which sets all pixel values above it to 
WHITE (level 255), and all those below to BLACK (level 0). The user then views 
the thresholded image and, if the result is not satisfactory, is given the option to 
enter a new threshold limit. Threshit.c utilizes a function named threshitf) which is 
physically located within the genfunc.c program file. This organization was due to 
the fact that other program files, such as feat-id. c , call the threshit() function. 

The key to the threshold process is a function call in the ITEX/PCp/us soft- 
ware package called THRESHOLD. The ITEX/PC function is also named THRESH- 
OLD and has the following format: 

void threshold (lowcut, highcut ); 

The lowcut and highcut integer values defined the lower and upper boundaries of the 
thresholded region. This threshold procedure could only be used on the currently 
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selected LUT [Ref. 12:p. 7-27]. The function in ITEX/PCp/us software used to set a 
selected LUT to a prescribed threshold is: 

void threshold(group, bank, highcut, lowcut)] 

where, 

group specifies the LUT: RED, GREEN, BLUE, or INPUT, 
bank specifies the output LUT used, in the range of 0 to 7, 
highcut specifies the upper bound for the threshold, and 
lowcut specifies the lower bound for the threshold. 

The threshit.c program utilizes all four LUT groups and the Oth 256-byte bank to 
specify a linear transform output. This specification differs significantly from the 
original ITEX/PC function in the amount of versatility granted the programmer. 
Any of the LUTs can be modified with this threshold command, not just the current 
LUT [Ref. 6:p. 7-35]. The lowcut value must never be specified as higher than highcut 
value, but these values may be made equivalent to display one intensity at level 255 
(WHITE). 

Once the user is satisfied with a particular image threshold, the image can 
be mapped to memory. This procedure involves the use of another ITEX/PCp/us 
function: 



int maplut (group, bank, x, y, dx, dy ); 

MAPLUT modifies actual memory values rather than just LUT values as THRESH- 
OLD does. The group and bank variables are identical to those used in the THRESH- 
OLD function, but now only the INPUT LUT is modified for mapping purposes in 
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threshit.c. The quantities x, y, dx, and dy refer to the horizontal and vertical coor- 
dinates of the 512x512 image [Ref. 6:p. 7-17]. As in the case of the THRESHOLD 
function, ITEX/PC permits the user to transform only that area specified by the 
current LUT. The variables group and bank are not defined in ITEX/PC MAPLUT 
function [Ref. 12:p. 7-14]. 

6. Program File feat-id. c 

It is mandatory that an image is thresholded prior to initiating the program 
feat-id. c (Appendix F). For this reason, the user is given the option within feat-id. c 
to threshold the image if thresholding was not completed previously. Conversion to 
a feature-labeled image requires an iterative process, whereby the image is evaluated 
row by row. Once this feature identification process is complete, the image can be 
further processed for specific feature data. 

The feat-id. c program necessitated few modifications for compatibility with 
ITEX/PCp/ws software. Again, the more extensive changes were performed on the 
general support functions. Feat-id. c used the ITEX/PCp/us functions WPIXEL 
and RPIXEL extensively, but these calls retained the same format as that used in 
ITEX/PC software [Refs. 6 and 12]. 

7. Program File sizeit.c 

Program sizeit.c (Appendix G) creates an output file size.dat which con- 
tains the dimensions of all image features in a format necessary for input into a sta- 
tistical analysis program. Sizeit.c relies upon pointer structures for dynamic memory 
allocation of processed feature data. 

A conversion factor, SCALE-FACTOR, figures prominently in this sizing 
algorithm. Kaeser determined this factor to be 3.7353 contingent on image calibration 
[Ref. 4], Upon subsequent runs of the sizeit.c program, however, a SCALE-FACTOR 



23 



of 10.0000 was necessary to achieve compatible results with Kaeser’s program [Ref. 
4:p. 17]. The constant term SCALE-FACTOR is defined in the header file thesis.h 
and can be easily modified if required. Based upon equipment configuration during 
hologram processing, the subject conversion factor is still questionable in terms of 
accuracy. For this reason, an additional function was devised to determine a con- 
version factor based upon user input. This function was entitled scale(), and its 
characteristics are discussed in the section relevant to genfunc.c. 

8. Program File speckle. c 

To calculate the value of speckle index (SI) after each filter iteration, 
speckle. c (Appendix H) is called to evaluate filtering effectiveness. This calculation 
is performed upon specified dimensions of the filtered image; integer values for im- 
age row and column are passed to this function from the filter routines. Values for 
speckle index typically range from a maximum value of 1.0, a rare occurrence, to 
a hypothetical value of 0 [Ref. l:p. 20]. No program modifications were necessary 
for ITEX/PCp/ws compatibility. Speckle index is obtained purely from statistical 
manipulations. 

9. Program File vir.arry.c 

Vir-arry.c (Appendix I) consists of a group of functions supporting setup 
and access of the virtual disk- based array stratagem employed primarily by the filter 
routines. The virtual array functions were adapted by Kaeser [Ref. 4] from Refer- 
ence 8. These arrays reside on disk and are accessed as though they are located in 
memory. This feature overcomes the memory limitation normally corresponding to 
array usage. Through the use of MACROS and pointer notation, the C programming 
language provides the means of managing virtual arrays [Ref. 8]. Thesis.h contains 
the obligatory MACROS needed for structure definition. All virtual array functions 
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mandated the use of long integers for indexing purposes [Ref. 13]. Since the dimen- 
sions of the actual image are defined in terms of integer values, a conversion interface 
problem developed in all the filter programs. This was further complicated by the 
fact that the ITEX/PCp/us WPIXEL function references datum points in terms of 
integer values. This function is essential for “writing” the filtered pixel values to 
memory from the virtual array routines [Ref. 6:p. 5-4]. With the use of the C 
language “(long)” command, the integer indices referencing pixel values within the 
image can be used simultaneously to reference values within the virtual arrays with 
no conversion difficulty [Ref. 15:p. 36]. 

10. Program File genfunc.c 

The program file genfunc.c (Appendix J) contains image processing general 
support functions. These functions maintain routine evolutions occurring within the 
main image processing programs. The following sections discuss functions that were 
either created or modified due to the advent of the ITEX/P Cplus libraries. 

a. Function startitQ 

This function makes provisions for initial Frame Grabber board setup. 
All the main function files with the exception of speckle. c and vir.arry.c access this 
function. Due to the introduction of the dual memory configuration, new ITEX/PCp/us 
calls were required in addition to what Kaeser utilized [Ref. 4]. The SETDIM rou- 
tine determines the configuration of the memory, either dual-store or single-store [Ref. 
6:p. 2-6]. SELECT_MEM determines which frame memory will be accessed while in 
the dual-store mode. Frame memory choices consist of MEM_A or MEM_B as de- 
fined in the itexpfg.h header file [Ref. 6:p. 2-19]. DISPLAY _MEM must be used in 
conjunction with the SELECT.MEM routine. It specifies which memory area will be 
displayed and gives the user the option of writing to one image area while displaying 
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the other [Ref. 6:p. 2-20]. For the current application, MEM_A is both the accessed 
and displayed frame memory. 

b. Function readit( ) 

The image is read from disk storage and placed on the video monitor 
via the readit( ) function. Any comment associated with the image when it was stored 
will be displayed at that time. The ITEX/PCp/us function READIM is the key to 
this procedure. The READIM call format is relatively unchanged from the ITEX/PC 
version. If an error occurs when a file is being read, the function will return an error 
code to identify why the operation failed. The itexpfg.h header file defines these error 
code values [Ref. 6:p. 6-5]. The standard ITEX/PC format did not make provisions 
for these values in its associated header file, itexpc.h ; therefore, file errors are referred 
to numerically [Ref. 12:p. 6-7]. Accordingly, the readit() function was modified to 
reflect the error nomenclature. 

c. Function saveitf) 

The saveitf) function allows the user to save the image to disk memory 
via the ITEX/PC plus function SAVIM. The image can be stored using compression; 
although compression may, in rare cases, cause the image to occupy more disk space 
than an uncompressed version. If compression is used, the file will be automatically 
decompressed when read from the disk. As in the case of the READIM function, 
SAVIM returns error code values if the save operation is unsuccessful. Subsequent 
modifications were performed to the saveit( ) code to make it compatible with the 
new documentation [Refs 6 and 12]. 

d. Function dev(stddev) 

Both program files Istat.c and 2sigma.c require the computation of the 
standard deviation of the image for their respective filtering algorithms. The function 
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dev(stddev) computes this value through an iterative process utilizing the RPIXEL 
function. Since the standard deviation is computed directly from the image without 
user input, there is a slight increase in filter program execution times associated with 
the implementation of this function. 

e. Function scale(factor) 

This function supports the sizeit.c program by allowing the user flex- 
ibility in the determination of a scale factor used in the sizing algorithm. The user 
is given the option of using the SCALE-FACTOR as delineated in the header file 
thesis.h , or of entering dimensions of a calibration object. If the user decides upon 
the latter, the length of the object must be entered in terms of pixels and inches. 
The function will then calculate the new SCALE-FACTOR and pass this value to the 
sizeit.c program. 

E. SUMMARY 

The modifications and improvements performed on the image processing pro- 
grams were necessary to make them compatible with the environment established by 
the installment of the PCVISIONp/us Frame Grabber board into the IBM PC/AT. 
These subject programs support all requirements for image acquisition, filtering, 
thresholding and particle identification. The following chapter examines the per- 
formance and verification of these programs at length. 
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IV. SOFTWARE PERFORMANCE ANALYSIS 



This chapter will discuss the pertinent details concerning software performance. 
The subjects of executable program size, timing analysis, and program verification 
will be examined at length in the following sections. Any departures from results 
obtained from the original C language programs will be appraised. 

A. EXECUTABLE PROGRAM SIZE 

Table 4.1 shows the comparison between executable program size of ITEX/PC 
and ITEX/PCp/us supported software. The size information is presented in terms of 
kilobytes (KB) and was acquired via the use of the DOS command, DIR [Ref. 4]. 
The ITEX/PC- based programs written by Kaeser [Ref. 4] have significantly smaller 
executable file sizes. This is attributed to the fact that during the LINK process, 
Kaeser linked main programs to required library files only. It was not necessary to 
link to the object modules of the genfunc.c, speckle. c, or the vir.arry.c files since all 
necessary functions were appended to the end of the main program files [Ref. 4]. In 
addition, the ITEX/PC function library differs from the ITEX/PCp/us version, and 
these differences will also impact executable file size. 

TABLE 4.1: Comparisons of Executable Program Size (KB) 



FILENAME 


ITEX/PC 


ITEX/PCp/us 


THRESHIT.EXE 


31,467 


63,101 


FEATJD.EXE 


32,509 


64,407 


SIZEIT.EXE 


46,281 


65,429 


2SIGMA.EXE 


49,4S3 


65,491 


LSTAT.EXE 


49,475 


65,499 


GEOFIL.EXE 


49,969 


66.317 



2S 



TABLE 4.2: Comparisons of Actual Program Size (KB) 



FILENAME 


ITEX/PC 


ITEX/PCp/us 


THRESHIT.C 


5,807 


1,336 


FEATJD.C 


10,130 


6,003 


SIZEIT.C 


11,011 


7,426 


2SIGMA.C 


14,623 


5,396 


LSTAT.C 


14,434 


5,396 


GEOFIL.C 


14,888 


7,916 



One stated objective of this thesis was to optimize the existing C language 
code. Although executable program size comparisons contradict this goal, it is shown 
in Table 4.2 that actual program sizes for the ITEX/PCp/us compatible programs are 
substantially reduced from the original versions. Two linker options, /F and /PAC, 
were initially employed to achieve more compact executable files. The use of these 
options should produce more rapid execution of files and shorter program load times 
[Ref. 16:p. 185]. The /E option had been utilized by Ivaeser in Reference 4, but 
this file packing option may not always give satisfactory results. It can occasionally 
increase file size, and this obviously reduces savings in disk space [Ref. 17:p. 268]. 

Upon linking program files with the /PAC and /F options, it was found that 
these options did not pose an improvement over the /E option. In fact, executable 
file sizes were increased. For example, the geofil.exe file was 73,132 KB when these 
options were used in the LINK process. The /E option produced an executable file 
size of 66,317 KB. All other executable files were affected similarly. Accordingly, the 
/E option was utilized for the LINK procedure. Appendix K contains the listings of 
the batch files used to compile and link the program files. 
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TABLE 4.3: Program Execution Times 



PROGRAM 


TIME(SEC) 


FEATURE IDENTIFICATION 


809.82 


FEATURE SIZING 


1780.03 


2SIGMA FILTER 


760.55 


LOCAL STATISTICS FILTER 


1324.47 


GEOMETRIC FILTER 


3403.13 



B. PROGRAM EXECUTION TIMES 

The jl7res4-img library image was used to test program execution times. This 
was the same image as that used by Kaeser in Reference 4, depicting an Air Force Res- 
olution Target (AFRT) obscured by speckle noise. The unfiltered image jl7resj.ing 
is depicted in Figure 4.1. Testing was accomplished by the use of a timer program, 
timer. c, in conjunction with a DOS input file to eliminate user intervention. Previous 
thesis work involved time testing of one quarter of the 512x512 image. MSFORT 
versions of the programs could effectively handle only one quadrant of the full screen 
image. Consequently, Kaeser conducted execution time comparisons on only one 
quarter of the image [Ref. 4]. Due to the banded nature of the ITEX/PCp/us frame 
memory, quadrant timing tests of the image would not produce results that could be 
compared with the original findings [Ref. 5]; therefore, timing analysis of the pro- 
grams was performed on the entire image. The resultant execution times are listed 
in Table 4.3. 

C. PROGRAM VERIFICATION 

To verify the results obtained from the ITEX/PCp/us-based programs, two im- 
ages were used from the existing image library. The jl7res4-img was again used, this 
time to determine the accuracy of the filter programs. Both featJd.c and sizeit.c were 
operationally tested with the image lOxwgrid.img derived from a standard calibration 
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Figure 4.1: Unfiltered jl7res4-img. 



image [Ref. 2]. Resultant data acquired from these images was compared to data in 
Reference 4. 

Speckle index value comparisons were used to determine the effectiveness of the 
filter programs. In each case, the respective filter was run through six iterations of 
the full screen image, jl7res4-img. Table 4.4 contains the tabulation of the results 
acquired, and the results from the original C programs are listed in Table 4.5. A file 
comparison was performed to ensure that the image used for testing was identical 
to the full screen image used previously. Comparison of the data obtained from the 
current iterations to data from Kaeser’s programs indicates slight variations in the 
speckle index values. The initial speckle index values for the two cases differ. One 
possible factor attributing to this difference is the inconsistency of integer conversions 
between the C language programs [Ref. 4]. The more significant measure of filter 
performance is the relative decrease in speckle index upon subsequent filter iterations. 
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TABLE 4.4: Speckle Index Reduction Results Obtained from Modified 
Programs 



NUM 


2 SIGMA 


L STAT. 


GEOMETRIC 


0 


0.293046 


0.293046 


0.293046 


1 


0.196282 


0.180061 


0.218982 


2 


0.144286 


0.125157 


0.174153 


3 


0.115754 


0.099063 


0.143760 


4 


0.098420 


0.084035 


0.121202 


5 


0.087068 


0.074063 


0.103756 


6 


0.079042 


0.066919 


0.089894 



TABLE 4.5: Speckle Index Reduction Results (Ref. 4) 



NUM 


2 SIGMA 


L STAT. 


GEOMETRIC 


0 


0.304005 


0.304005 


0.304005 


1 


0.153767 


0.194121 


0.224875 


2 


0.098967 


0.151951 


0.178222 


3 


0.074415 


0.129769 


0.146152 


4 


0.060615 


0.11593S 


0.122347 


5 


0.051691 


0.106372 


0.103983 


6 


0.045354 


0.099367 


0.089425 



Therefore, the dissimilarities of initial speckle index value are not detrimental to filter 
function, and all filter programs produce desirable reduction of speckle. Consequently, 
all evidence indicates that the ITEX/PCph/s-based filter programs function correctly. 
Figure 4.2 illustrates jllres^.img after six iterations of the geometric filter. 

The results obtained from the featJd.c and sizeit.c programs were identical to 
those obtained from Kaeser’s programs [Ref. 4]. The SCALE-FACTOR entry in the 
thesis. h header file had to be modified in order for the output files from the respective 
sizing programs to be congruent. Otherwise, the program files produced the same 
outputs as the original C language programs. 
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Figure 4.2: jl7res4.img Image After Six Iterations of the Geometric Filter. 
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V. CONCLUSIONS 



A. SUMMARY OF FINDINGS 

This thesis accomplished all major objectives as set forth in the introductory 
chapter. No major problems were anticipated with the integration of the YTEX/PCplus 
libraries and the COMPAQ 386/20 computer system. Due to a memory conflict cre- 
ated within the COMPAQ, the PCVISION plus had to be installed into the IBM 
PC/AT to complete C program modifications. Therefore, the bulk of this research 
was conducted on the PC/AT used in previous work. Although several of the pro- 
grams produced in Kaeser’s research [Ref. 4] were modified with no major difficulties, 
type definitions created major faults in the filter routines which had to be resolved. 
It was the goal of this thesis to improve program portability and optimization. For 
the most part, this goal was achieved. Executable program size did increase, but this 
was the result of program structure changes and library modifications introduced by 
the ITEX/Pcplus system upgrade. 

B. FUTURE WORK 

Further verification testing is necessary to sufficiently demonstrate the operation 
of the C language programs. A minimal amount of experimentation was accomplished 
to compare the acquired results with data obtained from the previous version of the 
programs. Further testing is mandatory for complete program validation and should 
be pursued in subsequent thesis research prior to porting the programs over to other 
computer systems. Some effort was expended on the application of the NDP C-3S6 
compiler to eliminate the need for virtual arrays for data processing. These efforts 
were largely unsuccessful due to the incompatibility of the ITEX/PC plus libraries 
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with the subject compiler system. This area of research could be further investigated 
for future thesis work. The following topics may also be considered for future work: 

• Utilization of Macintosh board for image processing. 

• Measurement of image particles utilizing a scanning electron microscope for 
sizing routine validation. 

• Further validation of image measurements using new calibration objects. 

• Integration of SUN workstations for image processing. 
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APPENDIX A: PROGRAM HEADER FILE 

thesis. h 



/* PURPOSE: 

To completely list all necessary include files, manifest constants, 
and MACRO definitions required by all the programs within this 
package. Any changes to constants in thesis. h will guarantee 
redefinition of those parameters throughout the main files. 



*/ 



/* Include files for use with ITEXPC programs 
*/ 



#include 

#include 

#include 

#include 

#include 

#include 

#include 

#include 



<stdio . h> 
<stdlib.h> 
<math.h> 
<malloc .h> 
<graph . h> 
<conio.h> 
"itexpf g . h" 
"stdtyp.h" 



/* Program constants used in the C program modules of the ITEXPC 

Optical System used for Rocket motor Hologram analysis. */ 

/* Intial ITEXPC Board Jumper Settings */ 

#define MEMBASE OxDOOOOL /* Itex board base memory start address */ 

#define REGBASE 0x300 /* Itex board base register start address */ 

#define MEMORY DUAL /* Itex board memory type */ 

/* ITEXPC initial AOI (area of interest) settings */ 



#def ine IXS 0 

#def ine IYS 0 

#def ine NROV 480 

#def ine NCOL 512 

#def ine COLOR 150 

#def ine XSIZE 512 

#def ine YSIZE 512 

#def ine DEPTH 8 



/* ItexPC LUT (Look Up Table) Variables */ 
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/* Threshold Limits */ 



#def ine 


LOWEST 


0 


/* Equates 


to Black 


for 


lowcut value */ 


#def ine 


HIGHEST 


255 


/* Equates 


to White 


for 


highcut value */ 


#def ine 


BLACK 


0 










#def ine 


WHITE 


255 











/* Filter array sizes and limits */ 

#def ine SNUM 25 

#def ine NUM 9 

#def ine HIGH 254 

#def ine LOW 0 

/* Sizing magnification factor used in sizing program */ 

#def ine SCALE_FACTOR 10.000000 /* Based on ( 317.5 / 85 ) */ 

/* Virtual array Header File Definitions */ 

#define header 7 



/* Virtual Array Control Block typedef */ 
typedef struct { 



FILE *f ile ; 
long size; 
int elsize; 
char * buffer; 
int buf_elsize; 
int buf_size; 
char *blank_rec; 



/* pointer to file control block */ 

/* number or array elements in file */ 

/* number or bytes in each element */ 

/* pointer to array buffer */ 

/* size of element in buffer including index */ 
/* number of elements in buffer */ 

/* pointer to initialization record */ 



> 

VACB 



/* used for extending file */ 

; /* Virtual Array control block type name */ 



/* Virtual Array Access Prototypes */ 

int init_v_array(char *f ilename , int rec_size , char filchar); 
VACB *open_v_array (char *f ilename , int buff er.size) ; 
void close_v_array(VACB *v_array) ; 
void *access_v_rec(VACB *v_array ,long index); 

/* Virtual Array Access Macros */ 

#define VREC(y) ((items *)access_v_rec(item_array , y) ) 
#define gdesc(y) VREC(y)->v_gdesc 
#def ine G(x,y) VREC(y)->v_gdesc [x] 

/* Virtual Array element structure typedef */ 
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typed©! struct { 

unsigned char v_gdesc [NCOL] ; 

> 

items ; 
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APPENDIX B: PROGRAM FILE: geofil.c 



/♦PURPOSE : Provides for filtering an image using the geometric filter 
algorithm. This program processes the image and provides 
for the operator to select the number of iterations. 

The area of filtering for the image is controlled by the 
operator by keyboard input of number of ROWS and COLS to 
process. Program will allow for a Maximum input of 480 rows 
and 512 cols with a minimum of 1 row and 1 column. 

The resulting filtered image pixel value is stored on 
disk in a virtual array, that is written to the image 
processing screen when all calculations are complete. 

Total calculation time for a 512x480 image is 3403.13 sec. 

DISK SPACE REQUIRED for virtual array storage is 250 KB minimum 
available on the disk at start of the routine. 



*/ 

#include "thesis. h" 

char f ilename[20] , comline [200] ; 

int srow, scol, times, row, col; 

main() 

{ 

int flag; 

extern int row, col, srow, scol, times; 

printf ("\n\tReady to Load IMAGE FILE?... Yes (y) / No (n) "); 
flag=getch(); 

if (flag == | | flag == 'y') 

{ 

startit () ; 
readit () ; 

> 

if ( flag == 'N ' || flag == 'n' ) return(0) ; 

geti("\n\nEnter number of rows (Max=480) to use in SPECKLE CALC : " ,&srow); 
geti("\nEnter number of columns (Max=512) to us in SPECKLE CALC : " ,&scol); 

geti("\n\n\nEnter number of filter iterations to do . . . " , Atimes) ; 
geti("\nEnter the number of IMAGE rows (MAX = 480) to FILTER : " ,&row); 
geti ("\nEnter the number of IMAGE cols (MAX = 512) to FILTER : ",&col); 
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switch(times) 

{ 

case 0: 

speckle() ; 

print! ( M \nSave image to Disk File?.. Yes (y) / Ho (n) "); 
flag=getch() ; 

if(flag == ’Y’ I | flag == 'y') saveitO; 

exit ( 1) ; 

break; 

case 1: 

break; 

default ; 

break; 

> 

do 

{ 

speckle() ; 
geof il() ; 

> 

while(times == 0); 

print! ( "\nSave image to Disk File?.. Yes (y) / No (n) "); 
f lag=getch( ) ; 

if (flag == >Y> I | flag == 'y') saveitO; 



/* geometric filter algorithm */ 



geof il() 



VACB 

extern int 

int 

int 



*item_array ; 

row, col, times; 

fl, f 2 , gl, g2 , g3 , x, y; 

a, b, c, d, pixel, calls; 



/* create a virtual array for the array of filtered 

pixel values the size of the image ( 480x512 ) */ 

init_v_array (" ITEMS .VAR" , sizeof (items) ,NULL) ; 

/* open the virtual array, reserve buffer space for 100 elements */ 



item_array = open. v_array ("ITEMS . VAR" , 100) ; 



/* start local statistic filter routine */ 



if ( row > NR0V ) row = NR0V; 
if ( col > NC0L ) col = NC0L; 
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/* Zero the initial array */ 



for( y = 1 ; y < row ; y++ ) 
for( x = 1 ; x < col ; x++ ) 

G( ( (long)x) , ( (long)y) ) = HULL; 

> 

> 

for(calls =• times; calls != 0; calls — ) 
c = 4; 

print! ("\n\n\tFilter Running — > */,d Runs alter this one ! . . . \n" , calls-1) ; 
for ( ; c ; ) 

{ 

switch(c) 

{ 

case 4 : 

a=l ; b=0; c=3 ; d=l ; 
print! ("\nc =• 3 "); 

tl (a,b, c ,d, item_array) ; 
break; 

case 3: 



a=0; 
b=l ; 
c=2 ; 

print! ("\nc = 3 "); 

tl(a,b, c,d,item_array) ; 
break; 

case 2: 



a=l ; 
b=l; 
c=l ; 

print! ("\nc = 2 *') ; 

tl(a,b,c,d,item_array) ; 
break ; 

case 1 : 

a=l ; 

b=-l; 

c=0; 

print! ( M \nc = 1 "); 
tl (a,b, c,d, it em_ array) ; 
printf("\nc = 0 END'*); 
speckleO ; 
break; 
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default : breaX; 

> 

> 

> 

/* close the virtual array */ 
close_v_array(item_ array) ; 



/* subroutine lor geofil ior steps 1 and 2 */ 

tl(a,b,c,d,it em_ array ) 
int a , b , c, d; 

VACB *item_array; 

{ 

int 11, 12, gl, g2, g3; 

int pixel, x, y; 

long xl, yl; 

extern int row, col, times; 

lor( y = 1 ; y < row-1 ; y++ ) 

{ 

lor( x = 1 ; x < col-1 ; x++ ) 

{ 

xl = ((long) x) ; 

yl = ((long) y); 

11 = rpixel( x,y ) ; 

12 = rpixel( x-a,y-b ); 

G(xl,yl)= max( 11, min( 12-1, 11+1 ) ); 

> 

> 

printl(" step 1"); 

lor( y = 1 ; y < row-1 ; y++ ) 

{ 

lor( x = 1 ; x < col-1 ; x++ ) 

xl = ((long) x); 
yl = ((long) y) ; 
gl = G( xl , yl) ; 

g2 = G(((long)(x-a)) , ((long) (y-b))) ; 

g3 = G(( (long) (x+a)) .((long) (y+b))) ; 

pixel = max( gl, min( g2, min( gl+1, g3+l ) ) ); 

wpixel(x,y, pixel); 

> 



printl(” step 2"); 
if ( d==l ) 

{ 

a= -a; 
b= -b; 
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d= 0; 

tl(a,b,c,d, item.array) ; 

> 

else il( d==0 ) 

{ 

d=l; 

t2(a,b,c,d,item_array) ; 
return; 

> 

> 



/* subroutine lor geolil lor steps 3 and 4 */ 
t2(a,b, c,d,item_array) 

int a, b, c, d; 

VACB *item_array ; 



{ 



int 

int 

long 

extern int 



fl, f 2 , gl, g2, g3; 

x, y, pixel; 

xi, yi; 

row, col, times; 



for( y = 1 ; y < row-1 ; y++ ) 

{ 

ior( x = 1 ; x < col-1 ; x++ ) 

{ 

xl = ((long) x) ; 
yl = ((long) y); 
fl = rpixel( x ,y ); 
f2 = rpixel( x-a,y-b ); 

G(xl,yl)= min( fl, max( f2-l, fl+1 ) ); 

> 

> 

printf(" step 3") ; 

for( y = 1 ; y < row-1 ; y++ ) 

{ 

for( x = 1 ; x < col-1 ; x++ ) 

{ 

xl = ((long) x) ; 
yl = ((long) y) ; 
gl = G( xl.yl); 

g2 = G(( (long) (x-a)), ((long) (y-b))); 

g3 = G( ( (long) (x+a) ) , ( (long) (y+b) ) ) ; 

pixel = max( gl, min( g2, min( gl+1, g3+l ) ) ); 

wpixel(x,y .pixel) ; 

> 

> 

printf(" step 4"); 
if ( d==l ) 
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a= -a; 
b= -b; 
d= 0; 

t2(a,b,c ,d, item_array) ; 

> 

else if ( d==0 ) 

{ 

d = 1; 
return; 

> 

> 
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APPENDIX C: PROGRAM FILE: lstat.c 



/♦PURPOSE : Provides for filtering an image using the local statistical 

algorithm. This program processes the image by using a local 
5x5 axray of pixels to calculate a statistical value for the 
central pixel of the local axray. 

The axea of filtering for the image is controlled by the 
operator by keyboard input of number of ROWS and COLS to 
process. Program will allow for a Maximum input of 480 rows 
and 512 cols with a minimum of 1 row and 1 column. Program 
also requires the image Std. Dev. as calculated by the genfunc.c 
Program. The resulting filtered image pixel value is stored on 
disk in a virtual axray, that is written to the image processing 
screen when all calculations axe complete. 

Total calculation time for a 512x480 image is 1324.47 sec. 

DISK SPACE REQUIRED for virtual axray storage is 250 FB 
minimum available on the disk at the staxt of the routine. 

*/ 

#include "thesis. h" 

char filename [20] , comline [200] ; 

int row, col, srow, scol, times; 

float devl; 

main() 

{ 

int flag; 

extern int row, col, times; 

printf ("\n\tReady to Load IMAGE FILE?... Yes (y) / No (n) "); 
flag=getch() ; 

if (flag == 'Y' | | flag == J yO 

{ 

startit () ; 
readit () ; 

> 

if (flag == 'N* || flag == J n') return(0) ; 

geti("\n\nEnter number of rows (Max=480) to use in SPECKLE CALC : ",&srow); 
geti("\nEnter number of columns (Max=512) to us in SPECKLE CALC : ",&scol); 

geti("\n\n\nEnter number of filter iterations to do . . . " ,&times) ; 
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geti("\nEnter the number o f IMAGE rows (MAX = 480) to FILTER : ",fcrow); 
geti("\nEnter the number of IMAGE cols (MAX = 512) to FILTER : ",ftcol); 
dev(ftdevl) ; 

printf ("\nThe standard deviation for this image is y,.2f\n" ,devl ) ; 



speckleQ ; 
lstat() ; 

> 

while(times == 0); 

printf ( "\nSave image to Disk File?.. Yes (y) / Ho (n) "); 
f lag=getch( ) ; 

if (flag == 'Y'll flag == 'y') saveitO; 

> 

/* */ 



Istat ( ) 



VACB 

long 

extern float 

extern int 

register int 

int 

float 

float 

int 



♦item.array ; 
xl.yl; 
devl ; 

row, col, times, srow, scol; 
m , n , nn j 

pixel ,a,b,x, y , calls ; 

suml , sum2, svar, lvar, tvar; 

Imean, lmean2, stddev, stddev2, k; 

ldata[SNUM] ; /♦ SNUM defined in Header file */ 



/* create a virtual array for the array of filtered 

pixel values the size of the image ( 480x512 ) */ 

init_v_array ( "ITEMS . VAR" , sizeof (items) , NULL) ; 

/* open the virtual array, reserve buffer space for 100 elements */ 
itera.array = open. v_array ("ITEMS . VAR" , 100) ; 

/* start local statistic filter routine */ 



if ( row > NR0W ) row = NR0W; 
if ( col > NC0L ) col = NC0L; 

stddev = devl / 255.0 ; /* Normalize deviation value */ 

stddev2 = stddev * stddev ; 

for(calls = times; calls != 0; calls — ) { 

printf ("\n\n\tFilter Running — > */,d Runs after this one ! . . . \n" ,calls-l) ; 
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for ( y = 2 ; y < row - 2 ; y++) 

{ 

lor ( x = 2 ; x < col - 2 ; x++) 

{ 

pixel = rpixel( x,y ); 
suml = 0; 
sum 2 = 0; 
nn = 0; 

for( n=y-2 ; n < y + 3 ; n++ ) 

{ 

lor( m = x- 2;m<x + 3; m++ ) 

ldata[nn] = rpixel( m,n ); 
suml += ldata[nn]; 
nn++ ; 

> 

> 

lmean = suml / SNUM ; 

lmean2 = lmean * lmean ; 

for( nn = 0 ; nn < SNUM ; nn++ ) 

{ 

svax - (ldata[nn] - lmean ) * ( ldata[nn] - lmean ); 
sum 2 += svar; 

> 

Ivar = sum 2 / SNUM ; 

tvar = fabs( ( (lvar + lmean2) / (stddev2 + 1) ) - lmean2 ) 
k = tvax / ( ( stddev2 * lmean2 ) + tvar ) ; 
xl = ((long) x) ; 
yl = ((long) y); 

G( xl,yl ) = (int) ( lmean + k * ( pixel - lmean ) ) ; 

} 

} 

/* write contents of the G array to image screen */ 

print! ( M \n\a\tWriting filtered image to screen...”); 

for( b = 2 ; b < row - 2 ; b++ ) 

{ 

for( a = 2 ; a < col - 2 ; a++ ) 

{ 

wpixel( a,b,G(((long)a) , ( (long)b)) ) ; 

} 

> 

speckle () ; 

/* close the virtual array */ 
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close_v_array(item_array) ; 



APPENDIX D: PROGRAM FILE: 2sigma.c 



/♦PURPOSE : Provides for filtering an image using a 2-sigma statistical 

algorithm. This program processes the image by using a local 
5x5 array of pixels to calculate a statistical value for the 
central pixel of the local array. The algorithm will eliminate 
pixels from the summation that are greater than 2-Sigma in 
value from that of the central pixel. 

The area of filtering for the image is controlled by the 
operator by keyboard input of number of ROWS and COLS to 
process. Program will allow for a Maximum input of 480 rows 
and 512 cols with a Minimum of 1 row and 1 column. Program 
also requires the image Standard deviation as obtained using 
genfunc.c program. The resulting filtered pixel value is 
stored on disk in a virtual array and written to screen later. 

Total calculation time for a 512x480 image is 760.55 sec 

DISK SPACE REQUIRED for virtual array storage is 250 KB 
minimum available on the disk at start of the routine. 



*/ 

#include "thesis. h" 



char filename [20] , comline [200] ; 

int row, col,t imes , srow, scol; 

float devl; 

main() 



int flag; 

extern int row, col, times; 



pr intf ( "\n\tReady to Load IMAGE FILE from disk?... Yes (y) / No (n) "); 
flag=getch() ; 

if (flag == 'Y' I I flag == 'y') 

{ 

startitO ; 
readit () ; 

> 

if(flag == , N I || flag == 'n') return(0) ; 

geti("\n\nEnter number of rows (Max=480) to use in SPECKLE CALC : ",&srow); 
geti("\nEnter number of columns (Max=512) to us in SPECKLE CALC : ",&scol); 

geti("\n\n\nEnter number of filter iterations to do. . . " ,&times) ; 
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geti("\nEnter the number o f IMAGE rows (MAX = 480) to FILTER : M ,*row); 
geti( M \nEnter the number of IMAGE cols (MAX = 512) to FILTER : ",*col); 
dev(ftdevl) ; 

print! ("\n\tThe standard deviation for this image is: 7.. 2f \n M ,devl ) ; 

do 

{ 

speckleQ ; 
sigmaO ; 

> 

while ( times == 0 ); 

print! ( M \nSave image to Disk File?.. Yes (y) / Ho (n) "); 4 

flag=getch() ; 

if (flag == 'Y'll flag == *y 1 ) saveitO; 



> 

/*******************************************************************/ 



sigraaO 

{ 

VACB 

extern int 
extern float 
long 

register int 

int 

int 

float 

int 



*item_array ; 

row, col t times, srow, scol; 
devl ; 
xl, yl; 
m, n, nn; 

pixel, suml , sum2, delta, high; 
low, a, b, x, y, calls; 
hvar, lvar, stddev; 

ldata[SNUM] ; /* SNUM defined in Header file */ 



/* create a virtual array for the array of filtered 

pixel values the size of the image ( 512x480 ) */ 

init_v_array( "ITEMS .VAR" , sizeof (items) , NULL) ; 

/* open the virtual array, reserve buffer space for 100 elements */ 
item_array = open_v_array ("ITEMS. VAR" , 100) ; 

/* start local statistic filter routine */ 



if ( row > NR0W ) row = NR0W; 
if ( col > NC0L ) col = NC0L; 

stddev = devl / 255.0 ; /* Normalize deviation value */ 

hvar = 1. +2. * stddev ; /* find high and low sigma */ 

lvar = 1. - 2. * stddev ; 

for (calls = times; calls != 0; calls — ) { 

printf ("\n\n\tFilter running — > 7d Runs after this one ! . . . \n" , calls-1) ; 
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for ( y = 2 ; y < row - 2 ; y++) 

{ 

for ( x = 2 ; x < col - 2 ; x++) 

{ 

pixel = rpixel( x,y ); 

suml = 0; 

sum 2 = 0; 

nn = 0; 

delta = 0; 

high = (int)(hvax * pixel); 
low = (int) (lvax * pixel); 

for( n=y-2;n<y+3; n++ ) 

{ 

for( m=x-2;m<x + 3; m++ ) 

{ 

xl = ((long) x); 
yl = ((long) y); 
ldata[nn] = rpixel( m,n ); 

if( ( low <= ldata[nn] ) ( ldata[nn] <= high ) ) 

l 

suml += ldata[nn] ; 
delta++; 

> 

nn++; 

> 

> 

if (delta <= 2) /* correct shot noise — 4 neighbor average */ 

{ 

sum2 = ( rpixel(x,y-l)+rpixel(x,y+l)+rpixel(x-l ,y)+rpixel(x+l ,y) ); 

G ( xl,yl ) = (pixel + sum2) / 5; 

continue; 

> 

G( xl,yl ) = suml / delta ; 

> 

> 

printf ( M \n\tWriting filtered image to screen..."); 

for( b = 2 ; b < row - 2 ; b++ ) 

{ 

for( a = 2 ; a < col - 2 ; a++ ) 

wpixel( a,b,G( ( (long) a) , ( (long)b) ) ); 

> 

speckleO ; 

/* close the virtual array */ 
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> 



close_v_air ay ( it em_ array ) ; 
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APPENDIX E: PROGRAM FILE: threshit.c 



/♦PURPOSE : This program thresholds the image on screen by taking 
an operator input value and forcing all image pixel 
values above the threshold value to the value of BLACK 
and all those below the value to the WHITE value. 

( Background == WHITE ; Feature == BLACK ) 

*/ 

#include "thesis. h" 

char filename [20] , comline[200] ; 

int srow, scol; 

mainQ 

{ 

int flag; 

printf ("\nThis program will threshold the input image desired."); 
printf ("\n\n\tReady to Load Image?... Yes (y) / No (n) "); 
f lag=getch() ; 

if (flag == 'YM I flag == 'y') 

{ 

startit() ; 
readit () ; 

> 

if (flag == 'N' | I flag == 'n') 

{ 

initializeQ ; 
return(O) ; 

> 

printf ("\n\n\tTHRESH0LD the image?... Yes (y) / No (n) "); 
flag = getch(); 

if (flag == 'YMI flag == 'y') threshitO; 



> 
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APPENDIX F: PROGRAM HEADER FILE 

feat Jd.c 



/♦PURPOSE : Labels and identifies each feature in an image on the ITEX 
system. Reads pixel-by-pixel and groups the objects by 
assigning a unique ID number to each feature or object so 
that they c an be processed by other routines. 

The ID routine requires a thresholded on image screen image 
or input of a previously saved thresholded image from disk or 
imput from disk and then thresholding. This program module 
includes an optional call to threshold if desired. 

( Background == WHITE ; Feature == BLACK ) 



*/ 

#include "thesis.h" 

char filename [20] , comline [200] ; 

int srow, scol; 

mainQ 

{ 

int flag; 

printf ( "\n\aYou MUST USE A THRESHOLDED IMAGE for this Program!!"); 
printf ("\n\n\tReady to Load image from disk?... Yes (y) / No (n) "); 
f lag=getch() ; 

if (flag == 'Y'll flag == J yO 

{ 

startit ( ) ; 
readit() ; 

> 

if (flag == 'N' || flag == 'n' ) return(O) ; 

printf ("\n\n\tNeed to THRESHOLD the image?... Yes (y) / No (n) ") ; 
flag = getch(); 

if (flag == 'YMI flag == ' y ' ) threshQ; 
f eat_id() ; 

printf ("\n\nSave image to Disk File?.. Yes (y) / No (n) "); 
f lag=getch( ) ; 

if (flag == 'YMI flag == 'y') saveitO; 

> 



54 



/* 


image feature 


identification algorithm */ 


feat 


.id() 






register int 


x, y; 




int 


a, xl, yl, x2 , y2, x3, y3, x4, y4; 




int 


fid = 0; 




int 


gid = 0; 




int 


temp_fid = 0; 




int 


temp_gid = 0; 




int 


maxfl, nl, n2, nla, nib, nlc, kk, kkl , nnl ; 



printf ("\n\n\tSTEP ONE in Progress..."); 

for( y = 1 ; y < NROW ; y++ ) 

for( x = 0 ; x < NCOL ; x++ ) 

if( rpixel(x,y) == WHITE ) continue; 
if ( x == 0 ) 

{ 

wpixel( x,y,fid ); 
continue ; 

> 

if( ( rpixel( x-l,y ) != WHITE ) ) 

{ 

wpixel( x,y ,rpixel(x-l,y) ); 
continue; 

} 

if ( rpixel(x,y-l) != WHITE ) 

wpixel( x,y ,rpixel( x,y-l ) ); 
continue ; 

} 

wpixel( x ,y ,f id ); 
a = x + 1 ; 

while (1) 

if ( rpixel( a,y ) == WHITE ) 

wpixel( x,y,fid ); 
step( ftfid.ftgid ); 
break; 

} 

if ( rpixel(a,y-l) != WHITE ) 

{ 

wpixel( x,y ,rpixel( a,y-l ) ); 
break; 

} 
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a += 1 ; 



> 

> /* end x - for */ 

} /* end y - lor */ 

fid += gid*HIGH; 

printf ("\n\n\tSTART STEP TWO... "); 

printf ("\n\nEnter the BEST estimate of the max feature length"); 
geti("\n\nlt * s Better to be too large than too small !!", toaxfl) ; 
printf ("\n\n\tRUNNING") ; 

for( y = 1 ; y < NROW ; y++ ) 

{ 

f or ( x = 1 ; x < NCOL ; x++ ) 

{ 

nl = rpixel( x,y ) ; 
n2 = rpixel( x,y-l ) ; 

if ( ( nl I | n2 == WHITE ) II (nl == n2 )) continue; 

x3 = x - maxfl; 
x4 = x + maxfl; 
y3 = y - maxfl; 
y4 = y + maxfl; 

if ( x3 < LOW) x3 = LOW; 
if ( x4 > NCOL ) x4 = NCOL; 

if ( y3 < LOW) y3 = LOW; 

if ( y4 > NROW ) y4 = NROW; 

for ( yl = y3 ; yl < y4 ; yl++ ) 

{ 

for ( xl = x3 ; xl < x4 ; xl++) 

{ 

nla = rpixel( xl.yl ); 

if ( nla == WHITE) continue; 

if ( nla == nl ) wpixel( xl.yl.n2 ); 

} 

fid -= 1; 

} 

> 

> 

printf ( "\n\n\tSTEP THREE in Progress"); 

for ( y = 1 ; y < NROW ; y++ ) 

{ 

f or ( x = 1 ; x < NCOL ; x++ ) 

{ 

kk = rpixel( x , y ) ; 

if ( kk == WHITE ) continue; 
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kkl = kk; 

nnl = temp_fid + temp_gid * HIGH; 
if( kkl < nnl ) continue; 
if ( (kkl - nnl) < -350 ) kkl = kk + 508; 
if ( (kkl - nnl) >= -350 ft& (kkl - nnl) < -150 ) 
kkl = kk + 254; 
if( kkl != nnl ) 

{ 

x3 = x - maxfl; 

x4 = x + maxfl; 

y4 = y + maxfl; 

if ( x3 < 0) x3 = LOW; 

if ( x4 > NCOL-1 ) x4 = NCOL -1; 

if ( y4 > HROW-1 ) y4 = NROW -1; 

f or ( yl = y ; yl < y4 ; yl++ ) 

{ 

f or ( xl = x3 ; xl < x4 ; xl++) 

{ 

nlc = rpixel( xl,yl ); 

if( ( nlc == WHITE ) I I ( nlc != kk ) )continue; 
wpixel( xl ,yl ,temp_f id ); 

} 

> 

step( &temp_f id ,Jttemp_gid ); 

> 

> 

> 

printf ( "\n\n\t\tFEATURE COUNT IS: */.d",fid); 

maplut ( INPUT, 0, IXS, IYS, NCOL, NROW ); 
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APPENDIX G: PROGRAM FILE: sizeit.c 



/♦PURPOSE : This routine uses an existing or saved image that has been 

thresholded and then ID'd with the ID Algorithm. The output 
of this program is a tabular output of the calculated Area, 
X-chord, and Y-chord that is suitable lor input to a 
Statistical Analysis Program. Output file is written to 
default drive as SIZE.DAT. 

( Background == WHITE ; Feature == BLACK ) 



*/ 

#include "thesis. h" 



char filename [20] , comline [200] ; 

int srow, scol; 



main() 

{ 

int flag; 

printf ("\n\aUSE IMAGE processed by the FEATURE ID program!!"); 



printf ("\n\tReady to Load ID'd image from disk?... Yes (y) / No (n) "); 
flag=getch(); 

if (flag == >Y> | | flag == 'y 1 ) 

{ 

startit () ; 
readit () ; 

> 

if ( flag == 'N' || flag == ) n I ) return(O) ; 

initializeO ; 
size_it () ; 



printf ("\n\tSave image to Disk File?... Yes (y) / No (n) "); 

f lag=getch() ; 

if (flag == >Y> | | flag == ' y 1 ) saveitQ; 



> 



/* image feature sizing algorithm */ 



size_it () 

{ 

register int nn, x, y, xa, ya, xb, yb; 

int j, ip, ipa, ipb, xl, yl, x2 , y2 , maxfl, num; 
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int 

int 

int 

int 

float 

float 

float 

long 

long 

long 

long 

FILE 



area, xchord, ychord, factor; 

fxmax, fymax, fxmin, fymin, fmaxarea, fminarea; 
cf_l, cf_2, mag; 



xmax, ymax, xmin, ymin, maxarea, minarea; 



xf lag=0 
yflag=0 
af lag=0 



♦aptr ; 
♦xptr; 
♦yptr; 
♦stop; 



♦featdata; 



geti("\nEnter the number of Features" , ftnum) ; 
getf ("\nEnter Magnification Factor" ,tonag) ; 
geti("\nEnter the Maximum Feature Length" ,ftmaxfl) ; 

/♦ Set initial test min/max values */ 

xmax = 0 ; 

ymax = 0; 

maxaxea = 0; 
xmin = 512; 

ymin = 512; 

minarea = 512; 

/* Dynamically allocate memory space for the number of 
features entered by the operator. 

Note: This eliminates the need to declare laxge arrays 
for these values and allows maximum number of features */ 

xptr = ( long *) calloc( num , sizeof (long) ) ; 

yptr = ( long *) calloc( num ,sizeof (long) ) ; 

aptr = ( long *) calloc( num , sizeof (long) ) ; 

/* Check for successful! memory allocation */ 



if( ! xptr || ! yptr II !aptr ) 

{ 

printf ("\n0ut of Memory!!! \n"); 

printf ("\nYou may have entered too many f eatures . \n" ) ; 
printf ("\nCorrect problem and try again."); 
exit (0) ; 



} 



/* Calculate Magnification Constants */ 



scale(fcfactor) ; 

cf_l = ( factor / mag ); 

cf _2 = ( cf_l ♦ cf_l ); 
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/* Above Conversion factor converts size to microns 

and is related to the input magnification for the 
hologram. The defined M SCALE_FACTOR ” constant 
value must be set for the proper value based on the 
equipment set up during actual image aquisition. 

The constant term SCALE.FACTOR is defined in the 
Header File THESIS. H 

*/ 

/* Begin Sizing Routine */ 

printf ("\n\n\t Star ting to SIZE!..."); 

nn = 0; /* Feature number counter */ 

f or( y = 0 ; y < NROV ; y++ ) 

{ 

if ( nn == num ) break; /* Quit when all features sized */ 

for( x = 0 ; x < HCOL ; x++ ) 

{ 

ip = rpixel(x , y) ; 

if ( (ip == WHITE) II (ip == BLACK) )continue; 

if( nn == num )break; /* Quit when all features sized */ 

xl = x - maxfl; /* Set up sizing Box */ 

yl = y; 

x2 = x + maxfl; 

y2 = y + maxfl; 

if ( xl < LOW) xl = LOW; 
if ( yl < LOW) yl = LOW; 
if ( x2 > HCOL ) x2 = NCOL; 

if( y2 > HROW ) y2 = NROW; 

f or ( ya = yl ; ya < y2 ; ya++ ) 

{ 

for( xa = xl ; xa < x2 ; xa++ ) 

{ 

ipa = rpixel(xa,ya) ; 
if ( ipa != ip ) continue; 
af lag++ ; 
xf lag++ ; 

> 

if ( xflag > *(xptr + nn ) ) *(xptr + nn) = xflag; 
xf lag = 0; 

> 

*(aptr + nn) = aflag; 
aflag = 0; 
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f or( xb = xl ; xb < x2 ; xb++) 

{ 

for( yb = yl ; yb < y2 ; yb++ ) 

{ 

ipb = rpixel(xb,yb) ; 

i f( ipb != ip ) continue; 

yflag++; 

wpixel(xb,yb, BLACK) ; 

> 

if ( yflag > *(yptr + nn ) ) *(yptr + nn) = yflag; 
yflag = 0; 

> 

/* Calculate Min/Max values */ 

xmax = max( xmax, *(xptr + nn) ); 

ymax = max( ymax, *(yptr + nn) ); 

max are a = max ( maxarea, *(aptr + nn) ) ; 
xmin = min( xmin, *(xptr + nn) ) ; 

ymin = min( ymin, *(yptr + nn) ) ; 

minarea = min( minarea, *(aptr + nn) ) ; 

nn++; /* increment counter */ 

printf ( H \n\n\tFeature V.d complete ... " ,nn) ; 

> 

} /*.... end initial x/y for loop ....*/ 

/* output data section * */ 

printf ("\n\a Sending output data to screen and FILE SIZE.DAT !"); 

/* wait(); */ 

if C ( featdata = f open("size .dat" , "w") ) == NULL ) 

{ 

printf ("CANNOT Open output file %s\n", *f eatdata) ; 
return(l) ; 

> 

printf (" ID NO AREA X-Vidth Y-Width \n M ); 

for( j = 0 ; j < num ; j++ ) 

{ 

f printf ( featdata, M, /.10d '/.10.3f y,10.3f 7.10. 3f\n", j+1, *(aptr + j) * cf_2, 
*(xptr + j) * cf_l, *(yptr + j) * cf_l ); 

printf ("\n%10d 7,10. 3f 7.10. 3f 7.10. 3f ", j+1, *(aptr + j) * cf_2, 

*(xptr + j) * cf_l, * (yptr + j) * cf _ 1 ); 

> 

/* This section can be used to put total information on bottom of data file 
if desired. File written now to be used with STATGRAPHICS */ 
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Max Y-chord= */.! 



Max Area= ’/•! " , 



/* !print! (!eatdata, "\n\n\nMax X-chord= '/.! 
xmax*c!_l ,ymax*c!_l ,maxarea*c!_2) ; 

!print! (!eatdata,"\nMin X-chord= */.! Min Y-chord= */,! Min Area= */•! ", 
xmin*c!_l ,ymin*c!_l ,minarea*c!_2) ; */ 

print! ( M \n\n\nMax X-chord= */.! Max Y-chord= */.! Max Area= */.! ", 
xmax*c!_l , ymax*c!_l ,maxarea*c!_2) ; 
print! ("\nMin X-chord= */,! Min Y-chord= */,! 
xmin*c!_l ,ymin*c!_l ,minarea*c!_2) ; 

/* Close open !ile and !ree allocated memory 

!close(!eatdata) ; 

!ree(xptr) ; 

!ree(yptr) ; 

!ree(aptr) ; 



Min Area= V.i ", 
*/ 
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APPENDIX H: PROGRAM FILE: speckle.c 



/♦PURPOSE : Provides for calculating the speckle index of an image. 

The resulting value is used as a measure to evaluate the 
sucess of filtering an image to reduce the speckle noise. 
This routine is called in all filter algorithms. 

— > OPERATOR must enter a Humber of ROWS and COLUMNS 
THE PROGRAM uses to calclate the speckle index. 

ROWS > Maximum = 480 Minimum = 1 

COLS > Maximum =512 Minimum = 1 

*/ 

#include "thesis. h" 
speckleO 

extern int srow, scol; 

int m, n, x, y, nn; 

long smax, smin; 

long sdata[NUM]; 

unsigned long deviation, ssum; 

float smean, slocal, stotal, tot; 
float spklindex; 

if ( srow >= NROW ) srow = NROW; 

if ( scol >= NCOL ) scol = NCOL; 

tot = (long)srow * scol; 

/* Commence calculation */ 

printf (”\n\n\tCalculating Speckle INDEX. . .”) ; 

for ( n = 1 ; n < srow-2 ; n++) 

{ 

for ( m = 1 ; m < scol-2 ; m++) 

{ 

smax = 0; 

smin = 260; 

ssum = 0; 

nn = 0; 
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for( y = ( n-1 ) ; y < (n+2) ; y++ ) 

{ 

!or( x = ( m-1 ) ; x < ( m+2 ) ; x++) 

{ 

sdata[nn] = rpixel( x,y ); 

ssum += (long)sdata[nn] ; 

i!( smax < sdata[nn] ) smax = sdata[nn] ; 

i!( smin > sdata[nn] ) smin = sdata[nn] ; 

nn++ ; 

> 

> 

deviation = smax - smin; 

if ( ssum == 0 ) smean =■ 1; 

smean = ssum / HUM ; 

slocal = (f loat)deviation / smean; 

stotal += slocal; 

> 

> 

spklindex = stotal / tot; 

print! ("\n\n\tThe Calculated '/,d by '/,d speckle index is '/,! ", srow , scol , spklindex) ; 

> 

/* END Speckle Function */ 
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APPENDIX I: PROGRAM FILE: vir_arry.c 



/♦PURPOSE : These functions provide for the setup of a disk drive virtual 
axray that can be indexed as if it were in the calling 
programs data storage area. The routines provide for disk 
file access when required to retrieve data elements from the 
axray. This allows the size of a declared array to be limited 
only by the amount of DISK SPACE available. See "thesis. h M 
header file for setup and MACRO routines that support these 
functions. Also program lstat.c fully implements the axray 
routines . 



♦/ 

#include "thesis. h" 

/♦ Virtual Array Access Routines ♦ / 

int init.v.axray (filename, rec.size .filchar) 

char ♦filename, filchar; 

int rec.size; 

{ 

long size; 

FILE ♦ f; 

f = fopen(f ilename,"wb") ; 
if (f ! = NULL) 

{ 

size = 0; 

fwrite(&size,4,l,f) ; 
fwrite(&rec_size,2, 1 , f ) ; 
fwrite(&f ilchar , 1, 1 ,f ) ; 
f close(f ) ; 
return(l) ; 

> 

else 

return(NULL) ; 

> 

VACB ♦open. v_array (filename, buff er_ size) 

char ♦filename; 

int buffer.size; 

{ 

VACB ♦v.array; 

char ♦buf.ptr; 

int i ; 

char filchar; 



/♦ write array size of zero ♦/ 
/♦ and array element size ♦/ 
/♦ and fill char ♦/ 
/♦ to file header ♦/ 
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/♦ allocate virtual array control block ♦/ 

v_array = (VACB ♦ ) malloc(sizeof (VACB)) ; 
if (v_ array == HULL) return(HULL) ; 

/♦ open virtual array file ♦/ 

v_array->f ile = f open(f ilename , ”r+b M ) ; 
if (v_array->f ile == HULL) 

f ree(v_array) ; 
return(HULL) ; 

>; 



/♦ get array size and element size for control block ♦ / 

f read(&v_array->size ,4, 1 , v_array->f ile) ; 
f read(Av_array“>elsize, 2 , 1 , v_array->f ile) ; 
f read(Af ilchar , 1 , 1 , v_array->f ile) ; 
v_array->buf _elsize = v_array->elsize + 4; 

/♦ allocate buffer ♦ / 



v_array->buf f er = (char ♦ ) malloc(v_array->buf _elsize * (buffer_size + 1)); 
if (v_array->buf f er == NULL) 

{ 

f close(v_array->f ile) ; 
free(v.array) ; 
return(NULL) ; 

>; 

v_array->buf_size = buffer_size; 

/♦ set up blank_rec using the fill character in array header ♦ / 

/♦ for initializing new array elements ♦/ 

buf_ptr = v_array->buf f er + v_array->buf _elsize * v_array->buf _size ; 

v_array->blank_rec = buf.ptr + 4; 

for (i =0; i < v_array->buf _elsize ; i++) 

♦buf_ptr++ = filchar; 

/♦ set record index negative for all buffer elements ♦/ 

buf_ptr = v_array->buff er ; 

for (i = 0; i < v_array->buf _size ; i++) 

{ 

♦((long ♦)buf_ptr) = -1L; 
buf_ptr += v_array->buf _elsize ; 

>; 

return(v_array) ; 
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> 



void 

VACB 

{ 



> 



void 

VACB 

long 

{ 



close_v_array(v_array) 
*v_array ; 



int i ; 

char *buf_ptr; 

long rec.index, file_offset; 

buf.ptr = v_array->buff er; 



/* Hush buffer */ 

for (i = 0; i < v_array->buf _size ; i++) 

{ 

/* check each element index */ 

/* if element present; write it to disk */ 

rec_index = *((long *)buf_ptr) ; 
if (rec_index >= 0) 

{ 

file_offset = header + rec_index * v_array->elsize ; 
f seek(v_array->f ile ,f ile_off set ,0) ; 

f write (buf_ptr + 4, v_array->elsize , 1 , v_array->f ile) ; 

>; 

buf_ptr += v_array->buf _elsize ; 

>; 

free(v_array->buff er) ; /* de-allocate buffer */ 

f close( v_array->f ile) ; /* close array file */ 

free(v_array) ; /* de-allocate VACB */ 



*access_v_rec( v_array , index) 
*v_array ; 

index ; 



char *buf_ptr; 

int buf_index; 

long rec.index, temp_ index; 



/* calculate buffer address of referenced element */ 
buf_index = index '/, v_array->buf _size ; 

buf_ptr = v_array->buff er + buf_index * v_array->buf _elsize ; 
rec_index = *(long *)buf_ptr; 
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/♦ if element present, return its buffer address ♦ / 
if (rec_index == index) return(buf _ptr + 4); 

/♦if element doesn't exist, extend the file */ 

if (index >= v_array->size) 

{ 

fseek(v_array->file,0,2) ; 

for (temp_index = v_array->size ; temp_index++ <= index; ) 

fwrite(v_array->blank_rec , v_array->elsize , 1, v_array->f ile) ; 
v_array->size = index + 1; 
f seek(v_array->f ile ,0,0) ; 

f write(&v_array->size , 4, 1, v_array->f ile) ; 

>; 

/♦ if buffer slot is occupied by another element, ♦/ 

/♦ save it to disk ♦ / 

if (rec_index >= 0) 

{ 

f seek(v_array->f ile , rec_index ♦ v_array->elsize + header, 0); 
f write(buf _ptr + 4, v_array->elsize , 1, v_array->f ile) ; 

>; 

/♦ read referenced element into buffer slot ♦ / 

f seek(v_array->f ile , index * v_array->elsize + header, 0); 
f read(buf _ptr + 4, v_array->elsize , 1, v_array->f ile) ; 

♦((long *)buf_ptr) = index; 

/♦ return address of element in buffer */ 

return(buf _ptr + 4) ; 



/♦ ♦/ 
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APPENDIX J: PROGRAM FILE: genfunc.c 



/♦PURPOSE : These general functions provide for routine evolutions 
that occur a number of times in the Image processing 
routines developed for analysis of speckle reduction 
algorithms. Some require support of other functions 
and axe not totally independant. 

*/ 

#include "thesis. h" 

/* initial ITEXPC Board setup */ 

startit () 

{ 

s ethdw (REGB ASE , MEMBASE , MEMORY ) ; 

setdim(XSIZE,YSIZE, DEPTH) ; 

aclear (IXS , IYS , NCOL , NROW , COLOR) ; 

fgon(); 

initial ize() ; 

cls() ; 

select_mem(MEM_A) ; 
display _mem(MEM_A) ; 

} 

/* Get integer keyboard input */ 

geti(name , iptr) 
char *name; 
int *iptr; 

{ 

printf(" %s M ,name); 
scanf (" '/,d M ,iptr); 

} 



/♦ Get floating point keyboard input */ 

getf ( r name, r iptr) 
char *mame; 
int *riptr; 

{ 

printf( M, /,s M ,rname); 
scanf (" */,f M ,riptr) ; 

} 



/* Pause routine used to prompt for operator intervention*/ 
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wait () 

{ 

print! ( M \n\n Press Return to continue"); 
f f lush(stdin) ; 
getcharO ; 

> 



/* returns minimum integer value from input array passed */ 

getmin(data , num) 

int *data; 
int num; 

{ 

int t, min; 

for (min= data[0] , t=l ; t<num ; t+ + ) 

if ( data[t] < min ) min = data[t] ; 

return min; 

> 

/* returns maximum integer value from input array passed ♦/ 

getmax (data, num) 

int *data; 
int num; 

{ 

int t , max ; 

for (max= data[0], t=l ; t<num ; t++) 

if ( data[t] > max ) max = data[t] ; 

return max; 

} 

/* DOS System Call to clear Display Screen */ 

cls() 

{ 

_clearscreen(_GCLEARSCREEN) ; 

> 

/* routine to read an image from disk */ 
readit () 

/♦program uses the following externals: filename , comline*/ 

{ 

extern char filename [20] , comline [200] ; 
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int 

int 



errval; 

ch; 



while(l) 

{ 

printf ("\n\nENTER IMAGE FILENAME as (DEV : name . IMG) ->") ; 

fflush(stdin) ; 

scanf ("*/,s" .filename) ; 

errval = readim(IXS, IYS ,NCOL,NROW, filename , comline) ; 
if (errval == 0) { 

printf ("\nFILE-> */.s \n\nCOMMENT: */,s" .filename , comline) ; 
break; 

> 

if (errval != 0) { 

switch (errval) 

case FILE_ERROR: 

printf ("Error opening file\n"); 
break; 

case F0RMAT_ERR0R: 

printf ("Unknown file format\n"); 
break; 

case READ_ERR0R: 

printf ("Error Reading file\n"); 
break; 
default : 

printf ("? Unknown Error ?’/.d\n" , errval) ; 
break ; 

> 

printf ("Try Again?? . . .Yes(y) / No(n)\n M ); 
ch = getch() ; 

if (ch == 'Y* || ch == ’y } ) continue; 
if(ch == 9 N , || ch == 'n') exit(0); 

> 

> 

} 



/* routine to save an image to disk */ 
saveit () 

/♦program uses the following externals: *f ilename , *comline*/ 

{ 

extern char f ilename [20] , comline [200] ; 

int format, errval; 

int ch ; 

initializeO ; 
while (1) 

{ 
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print! ("\n\nENTER FILENAME . . . (DEV : name . IMG)->") ; 

ff lush(stdin) ; 

scam! ("'/*s" ,f ilename) ; 

geti("\n Enter File Compression Value ? . . . (0/1) " ,*format ) ; 
print! ( M \n\n Enter Image Comment ..to continue.. (Max 200 CHAR)..."); 
f f lush(stdin) ; 
gets(comline) ; 

if (format == l)printf ("\n\nStoring Image Using Compression — Please wait 
if(format == 0)pr intf ( M \nStoring Image — Please wait!!"); 

errval = saveim(IXS, IYS,NCOL,NROV, format , filename ,comline) ; 

if (errval == 0) { 

print! ("\nlmage save completed Sat"); 
break; 

> 

if(errval != 0) { 

printf("\n\a Error saving file!!"); 
if (errval == ALL0CATI0N.ERR0R) { 

print! ("\nlnsuf f icient Disk Space"); 

> 

if (errval == WRITE.ERROR) { 

print! ("\nError writing file or values"); 

> 

print! ("\n\tTry Again?? . . .Yes(y) / No(n) M ); 
ch = getch() ; 

if(ch == || ch == J yO continue; 

if(ch == 'N 1 || ch == ’n’) break; 

> 

> 

> 

/* counter step routine */ 

step(u , v) 
int *u,*v; 

{ 

*u += 1 ; 

if (*u >= HIGH) { 

*u = 1 ; 

*v += 1; 

> 

> 

/* threshold routine for feat_id.c*/ 

thresh( ) 

{ 

int ans , c , option; 
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unsigned int val; 



while( 1) 

{ 

initluts () ; 
initializeQ ; 

get i ( ”\n\n\t ENTER HEW THRESHOLD LEVEL (0-255) :”,*v al) ; 
setlut ( INPUT, 0) ; 
threshold(INPUT,0, HIGHEST, val) ; 
setlut(RED,0) ; 

threshold(RED , 0 , HIGHEST , val ) ; 
setlut (GREEN, 0) ; 
threshold (GREEN , 0 , HIGHEST , val) ; 
setlut(BLUE,0) ; 

threshold(BLUE,0, HIGHEST, val); 
setlut (INPUT, 0); 

printf ( M \n\nSATISFIED WITH THRESHOLD AT 7.d? . . . Yes (y ) / No(n) . . . " , val) 
ans = getchO ; 

if (ans == ' Y ; II ans == *y ; ) break; 

> 

maplut (INPUT , 0 , IXS , IYS , NCOL , NROW) ; 
initializeQ ; 



/* threshold routine for threshit.c */ 

threshit () 

{ 

int ans , c , option; 
int *val; 

while(l) 

{ 

initluts () ; 
initializeO ; 

geti( M \n\nENTER NEW THRESHOLD LEVEL (0-255) : M ,&val); 

printf ("’/•d\n" ,val) ; 

setlut (INPUT , 0) ; 

threshold( INPUT , 0 , HIGHEST , val) ; 

setlut (RED, 0) ; 

threshold(RED,0, HIGHEST, val) ; 
setlut (GREEN, 0) ; 
threshold (GREEN , 0 , HIGHEST , val) ; 
setlut(BLUE.O) ; 

threshold(BLUE,0, HIGHEST, val); 
setlut (INPUT , 0) ; 

printf ( M \n\nSATISFIED WITH THRESHOLD AT */.d? . . . Yes (y) / No(n) . . . M , val) 
ans = getchQ ; 
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if (ans == 1 Y ’ || ans == 'y') break; 



> 

printf (”\n\t\t 0: SAVE THE THRESHOLD IMAGE?...”); 

printf (”\n\t\t 1: QUIT (Leave thresholded image for further processing)...”) 
printf (”\n\t\t 2: RESTORE system to original input image...”); 
geti(”\n\t Select option by NUMBER: ”,*option); 
switch(option) 

{ 

case 0: 

maplut ( INPUT , 0 , IXS , IYS , NCOL , NROW) ; 
saveit () ; 
break; 
case 1: 

maplut ( INPUT , 0 , IXS , IYS , NCOL , NROW) ; 
break; 
case 2: 

initializeO ; /*initialize without mapping operation*/ 

break; 
default : 
break; 

> 

> 

/* function to compute standard deviation of image */ 

dev(stddev) 

float *stddev; 

{ 

float var, sum, sqsum, sqvalue, dim, value; 

long x, y; 

sum = sqsum = 0; 

dim = (float) NROW * NCOL; 

f or( y = 0 ; y <= NR0W-1 ; y++) 

{ 

f or ( x = 0 ; x <= NC0L-1 ; x++) 

{ 

value = rpixel(x,y); 
sqvalue = value * value; 
sum += value/dim; 
sqsum += sqvalue/dim; 

> 

> 

var = sqsum - (sum* sum) ; 

*stddev = sqrt(var) ; 



/* function to compute scale factor for sizing routine */ 
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scale(factor) 

float 



♦factor; 



int inch, pixel, flag; 

float dinch, dpixel; /♦ default length values ♦/ 

/♦ dinch and dpixel axe to be modified accordingly ♦/ 
dinch = 1.00; 

dpixel = 50.00; 

printf ( M \n\tUtilizing standard SCALE_FACT0R? . . .Yes (y) / No (n)"); 
flag = getch(); 

if (flag == I I flag == >y') 

{ 

♦factor = ( (f loat)SCALE_FACT0R) ; 

} 

if (flag == , N > | | flag == 'nO 

{ 

printf ("\n\tPress 'O' if using default values for length in inches or pixels") 
geti("\n\tEnter the length of object in inches (default) ,&inch) ; 
geti("\n\tEnter the length of object in pixels(def ault) : " ,&pixel) ; 

if (inch == 0) 

{ 

inch = dinch; 

> 

if (pixel == 0) 

{ 

pixel = dpixel; 

} 

♦factor = ((float)inch) / ( (f loat)pixel) ; 

} 

printf ("\n\tSCALE_FACT0R is '/,f \n" , *f actor) ; 

> 
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APPENDIX K: C COMPILER OPERATIONS 



The following batch files list the C compiler directives used to compile and link 
all C programs. The C compiler files ACCOMP.BAT and DACCOMP.BAT are 
listed below. ACCOMP.BAT contains the switches necessary to utilize CODEVIEW 
debugger, and ACCOMP2.BAT is for final program optimization. The F drive listed 
in the PATH statement of these batch files is due to the fact that the programs were 
stored on a Bernoulli disk. 



ACCOMP.BAT: 

path c : \; c : \dos ; c: \ws jet ; f:\kedit.*; c:\mous el ;c:\xtree; c:\lsrctrl; . . . 
i :\;f :\msc\bin;f :\msc\tmp;f :\msc\src;f :\msc\wrk;f :\msc\lib; . . . 
f : \msc\include 
set include=\msc\include 
set lib=\msc\lib ; \pcplus\itex\lib 
set tmp=\msc\tmp 

cl /c /Fs /Zi /GO /AL /Fm B /.l.c 

ACC0MP2.BAT 

path c :\; c : \dos ; c:\wsjet ; f :\kedit . * ; c: \mousel ; c : \xtree; c: \lsrctrl ; . . . 
f : \;f : \msc\bin;f :\msc\tmp;f :\msc\src;f :\msc\wrk;f :\msc\lib; . . . 
f : \msc\include 
set include=\msc\include 
set lib=\msc\lib;\pcplus\itex\lib 
set tmp=\msc\tmp 

cl /c /Fs /Ze /GO /AL V.l.c 



The following batch files, LTHESIS.BAT and DTHESIS.BAT, link the programs 
to the libraries and supporting function files. The speckle. c should not be included 
within these files when linking the programs sizeit.c and feat.id.c due to an undefined 
external problem associated with the external variables SROW and SCOL. LTHE- 
SIS.BAT is used for debugging applications, and DTHESIS.BAT is for final program 
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optimization. It is also unnecessary to link feaLid.c and sizeit.c to the virtual array 
functions contained within vir.arry.c. 



LTHESIS.BAT: 

IF NOT ERRORLEVEL 1 LINK /NOD /CO /E ’/.l genfunc speckle vir.arry . . . 
, , , ITEXPCML . LIB LLIBCE . LIB ; 



DTHESIS.BAT: 

IF NOT ERRORLEVEL 1 LINK /NOD /E */.l genfunc speckle vir.arry... 

, , , ITEXPCML . LIB LLIBCE . LIB ; 

To create the menu format for program display, a C program named hologram. c 
was generated. To link all image processing programs to this module, all program 
“main()” statements were modified so that hologram. c was made the main program. 
All image processing programs were then linked to hologram. c after compilation. The 
ACCOMP2.BAT file was used for compilation, and the PTHESIS.BAT file was used 
for subsequent module linking. Due to the length of the LINK statement, the user 
must enter the libraries ITEXPCML and LLIBCE manually into the keyboard each 
time the LINK process is accomplished. 

PTHESIS.BAT: 

IF NOT ERRORLEVEL 1 LINK /NOD /E */.l table speckle vir.arry... 

genfunc threshit sizeit feat.id geofil.. 
lstat 2sigma 
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